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Abstract. The Polya number of a classical random walk on a regular lattice is known 
to depend solely on the dimension of the lattice. For one and two dimensions it equals 
one, meaning unit probability to return to the origin. This result is extremely sensitive 
to the directional symmetry, any deviation from the equal probability to travel in each 
direction results in a change of the character of the walk from recurrent to transient. 
Applying our definition of the Polya number to quantum walks on a line we show that 
the recurrence character of quantum walks is more stable against bias. We determine 
the range of parameters for which biased quantum walks remain recurrent. We find 
that there exist genuine biased quantum walks which are recurrent. 
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1. Introduction 

Random walks are a popular topic in physics i2j. The popularity stems from several 
sources. First, random walks are rather simple in their formulation yet powerful in their 
application and allow to pinpoint the essential physics involved in the studied processes. 
Next, the random walks are one of the tools which allow to connect the microdynamics 
with the macrobehaviour of large systems. Finally, random walks are quite flexible and 
popular also outside physics to describe various phenomena. Hence it comes not as a 
surprise that the first random walks have not been formulated within physics but to 
describe the alternation of share prices on the stock exchange or the spreading of insects 
in a forest [3l H] . 

The study of random walks obtained a new stimulus when they were combined 
with quantum mechanics [5l [6]. Here the walker is thought to be a non-classical object 
enriched with wave attributes. The novel features of quantum walks have been shown to 
be not only of theoretical interest but to also have practical implications, especially for 
quantum algorithms [71 [U [9l [10] . An important concept is the hitting time [TTl [T2l [T3[ [T^ 
which helps to point out the fundamental difference between classical and quantum 
walks allowing for algorithmic speed-up. One of the simplest non-trivial examples for 
a quantum walk is the one on a line [15] which is closely related to the so called 
optical Galton board [I6|. Various aspects of one dimensional quantum walks have 
been analyzed [TTl [HI [191 [SHI [211 [22]. Additional interesting effects, e.g. localisation, 
arise when one considers multi-state quantum walks [231 [211 [23 [26] . 

One of the characteristics of the random walk on an infinite lattice is expressed by 
the probability of the walker to return to its starting position, called the Polya number 
|27] . If the Polya number equals one the walk is called recurrent, otherwise there is a 
non-zero probability that the walker never returns to its starting position. Such walks 
are called transient. The recurrence nature of the random walk is determined by the 
asymptotic behaviour of the probability at the origin [2H]- One finds that a random 
walk is transient if the probability at the origin decays faster than t~^. The recurrent 
behaviour has been studied in great detail for classical random walks in dependence on 
the dimension and the topology of the lattice [291 130] • 



Recently, we have extended the concept of Polya number to quantum walks [3T] . In 
our definition we proposed a particular measurement scheme to minimize disturbance: 
each measurement in a series is carried out on a different member of an ensemble of 
equally prepared quantum systems. We have shown that the recurrence nature of 
the quantum walk, according to the above definition, is determined by the asymptotic 
behaviour of the probability at the origin in a similar way as in classical random walks. 
However, due to interference the asymptotics of the probability at the origin does not 
depend solely on the dimension of the lattice, but also on the coin operator and the initial 
coin state. Hence, one can find strikingly different recurrence behaviour for quantum 
walks compared to their classical counterpart [32]. Note that recurrence is meant here 
as the return to the origin which can be considered as a fractional recurrence from the 
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point of view of the whole quantum state [221 El] So far we have considered balanced 
walks, i.e. ones where there is no preference in direction for the walker and the step 
lengths are equal. For a large class of quantum walks this assumption does not hold 
and we wish to study the implications of unbalanced coins and unequal step lengths for 
the recurrence properties. 

In the present paper we study biased quantum walks on the line and compare 



their properties with their classical counterparts. As we briefly review in Appendix A 



recurrence of classical random walks is a consequence of the walk's symmetry. They are 
recurrent if and only if the mean value of the position of the particle vanishes. This 
is due to the fact that the spreading of the probability distribution of the position is 
diffusive while the mean value of the position propagates with a constant velocity. In 
contrast, for quantum walks both the spreading of the probability distribution and the 
propagation of the mean value are ballistic. We show that this allows for maintaining 
recurrence even when the symmetry is broken. 

Our paper is organized as follows: In Section[2]we describe the biased quantum walk 
on a line. In Section [3] we solve the time evolution equations with the help of the Fourier 
transformation. We find that the probability amplitudes can be expressed in terms of 
integrals where time enters only in the rapidly oscillating phase factor. This fact allows 
a straightforward asymptotic analysis of the probability amplitudes by means of the 
method of stationary phase. We perform this analysis in Section HI Since the recurrence 
of the quantum walk is determined by the asymptotics of the probability at the origin we 
find a condition under which the biased quantum walk on a line is recurrent. In Section [5] 
we analyze the recurrence of biased quantum walks from a different perspective. We find 
that the recurrence is related to the velocities of the peaks of the probability distribution 
generated by the quantum walk. The explicit form of the velocities leads us to the same 
condition derived in Section HI Finally, in Section [6] we analyze the formula for the 
mean value of the position of the particle derived in Appendix B in dependence of the 
parameters of the walk and the initial state. We find that there exist genuine biased 
quantum walks which are recurrent. Conclusions and outlook are left for Section [71 



2. Description of the walk 

Let us consider biased quantum walks on a line where the particle has two possibilities 
— jump to the right or to the left. Without loss of generality we restrict ourselves to 
biased quantum walks where the jump to the right is of the length r and the jump to 
the left has a unit size. We depict the biased quantum walk schematically in Figure [H 
The Hilbert space of the particle has the form of the tensor product 

n = np®nc (i) 

of the position space 

Hp = f (Z) = Span{|m)| m e Z}, (2) 
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Figure 1. Schematics of the biased quantum walk on a Une. If the coin is in the state 
\R) the particle moves to the right to a point at distance r. With the coin state \L) 
the particle makes a unit length step to the left. Before the step itself the coin state is 
rotated according to the coin operator C{p). 



and the two dimensional coin space 

^c = C2 = Span{|i?),|L)}. (3) 
A single step of the quantum walk is given by the propagator 

U = SiIp®C). (4) 

Here Ip denotes the unit operator acting on the position space Hp. The displacement 
operator S has the form 

+00 +00 

S = \m + r){m\®\R){R\+ ^ |m - 1) (m| ® |L) (L|. (5) 

m=— oo m=— oo 

The coin flip C is in general an arbitrary unitary operator acting on the coin space Tic 
and is applied on the coin state before the displacement 5* itself. However, as has been 
discussed in [17] the probability distribution is not affected by the complex phases of 
the coin operator. Hence, it is sufficient to consider the one-parameter family of coins 



From now on we restrict ourselves to this family of coins. The value of p = 1/2 
corresponds to the well known case of the Hadamard walk. 
We write the initial state of the particle in the form 

+00 L 

1^(0))^ J2 5^^.(m,0)|m)®K). (7) 

m=— oo i=R 

The state of the walker after t steps is given by successive application of the time 
evolution operator given by Eq. (jlj) on the initial state 

+00 L 

\^{t))^U'\m)= Yl Y.'^^M\m)^\^). (8) 

m=— oo i=R 

The state of the particle is fully determined by the set of two-component vectors 

tp{m, t) = (V^ij(m, t), i)L{m, t) f. (9) 
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Here ipR(^L) ij^i t) is the probability amplitude to find the particle at position m after t 
steps with the coin state \R{L)). The probability distribution generated by the quantum 
walk is given by 

P{m,t) = \{m,Rm))\' + \{m,L\i;m' 

(10) 



3. Time evolution of the walk 



To obtain explicit and closed form expressions for the time dependent state vector we 
rewrite the time evolution equation ([8]) for the state vector \ip(t)) into a set of difference 
equations 

ip{m, t) = C^{p)ip{m — r, t — 1) + 

+ C„(p)^(m + l,t-l) (11) 

for the probability amplitude vectors il){m,t). The form of the matrices C±{p) follows 
from the matrix C{p) 




The time evolution equations (ITT!) are greatly simplified with the help of the Fourier 
transformation 

+ 00 

ij{Kt)^ J2 ^(m,t)e*'"^ (13) 

m=— oo 

where the momentum is a continuous parameter ranging from — vr to n. The new 
function ip{k^t) is square integrable on a unit circle. 

The time evolution in the Fourier picture turns into a single difference equation 

'4){k,t) = U{k)^{k,t-l), (14) 

where the propagator has the form 

i/pe*'''" a/1 - pe*'''' 
a/1 - pe"*^ -^fpe 

The solution of (|T4l) is straightforward. We find 

'4,{k,t) = U\k)^{k,Q), (16) 

where ^/'(/c,0) is the Fourier transformation of the initial state. We restrict ourselves 
to the situation where the particle is initially localized at the origin as dictated by the 
nature of the problem we wish to study. As follows from (fT3l) the Fourier transformation 
^(/c, 0) of such an initial condition is equal to the initial state of the coin 



mk)^\ 1. (15) 
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which we denote by ip. Since ip can be an arbitrary normahzed complex two- component 
vector we parameterize it by two parameters a G [0, 1] and cp G [0, 27r) in the form 



'a 



To evaluate the powers of the propagator U{k) it is convenient to diagonalize it. Since 
the propagator is unitary its eigenvalues have the form e*'^^'^ where the phases read 

Ui{k) = — - — k + arcsin ^^psin ^ ^ fc^ ^ , 

r-1 / f r + 1 \\ 

^2{k) = — - — k — TT — arcsin ( y^sin ( — - — k j 1 (19) 

We denote the corresponding eigenvectors by vi^2{k)- We give their explicit form in the 



Appendix B With this notation we write the solution of the time evolution equation 
in the Fourier picture in the form 

2 

4^{k,t) = J2e''''^'^'iv,{k),^)v,{k). (20) 
i=i 

Here (, ) means scalar product in the coin space. Finally, we obtain the solution in 
position representation by performing the inverse Fourier transformation 

2 



Hk ~ rlh 

— ^{k,t) = Y. — e*('^^('^)*-"'=) {v,{k),^)vj{k). 

(21) 



4. Asymptotics of the quantum walk and recurrence 

To determine the recurrence nature of the biased quantum walk we have to analyze the 
asymptotic behaviour of the probability at the origin [31] . Exploiting fl2T|) the amplitude 
at the origin reads 

Hk 

mt) = Y.J ^e^^^^'^' {v,{k),^)v,{k), (22) 

which allows us to find the asymptotics of the probability at the origin with the help 
of the method of stationary phase [35] . The important contributions to the integrals in 
(122]) arise from the stationary points of the phases 0191) . We find that the derivatives of 
the phases uji^2{k) are 



r 

uj'^ik) = - 



2 ^4 + 2p [cos(A;(r + 1)) - 1] ' 
1 ^(r + 1) cos (A;^) 



v/4 + 2p [cos(A;(r + 1)) - 1]' 

(23) 
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Using the method of stationary phase we find that the amphtude will decay slowly - 
like t~^, if at least one of the phases has a vanishing derivative inside the integration 
domain. Solving the equations uj[ 2{k) = we find that the possible saddle points are 



2 / (l-p)(r-l)2\ ^ ^ 

ko = ± arccos ± W ^ — . 24 

r + 1 \ V P''^ I 

The saddle points are real valued provided the argument of the arcus-cosine in f l21j) is 
less or equal to unity 

^ '-^ < 1. 25 

Apr ~ 

This inequality leads us to the condition for the biased quantum walk on a line to be 
recurrent 

PH(r)>(^)'. (26) 
We illustrate this result in Figure [2] for a particular choice of the walk parameter r = 3. 



ex. 0.5 - 

L 




-TT -f f 

k 

Figure 2. The existence of stationary points of the phases wi_2(fc) in dependence on 
the parameter p and a fixed step length r. We plot the imphcit functions uj'i 2(fc) = 
for r = 3. The plot indicates that for p < p_r(3) = j the phases u}i^2{k) do not have 
any saddle points. Consequently, the probability amplitude at the origin decays fast 
and such biased quantum walk on a line is transient. For p > pr(S) the saddle points 
exist and the quantum walk is recurrent. 



Our simple result proves that there is an intimate nontrivial link between the length 
of the step of the walk and the bias of the coin. The parameter of the coin p has to be 
at least equal to a factor determined by the size of the step to the right r for the walk 
to be recurrent. We note that the recurrence nature of the biased quantum walk on a 
line is determined only by the parameters of the walk itself, i.e. the coin and the step, 
not by the initial conditions. The parameters of the initial state a and (p have no effect 
on the rate of decay of the probability at the origin. 
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5. Recurrence of a quantum walk and the velocities of the peaks 

We can determine the recurrence nature of the biased quantum walk on a line from a 
different point of view. This approach is based on the following observation. The well 
known shape of the probability distribution generated by the quantum walk consists of 
two counter-propagating peaks. In between the two dominant peaks the probability is 
roughly independent of m and decays hke t~^. On the other hand, outside the decay 
is exponential as we depart from the peaks. As it has been found in [15] the position 
of the peaks varies linearly with the number of steps. Hence, the peaks propagate with 
constant velocities, say vl and vr. For the biased quantum walk to be recurrent the 
origin of the walk has to remain in between the two peaks for all times. In other words, 
the biased quantum walk on a line is recurrent if and only if the velocity of the left peak 
is negative and the velocity of the right peak is positive. 

The velocities of the left and right peaks are easily determined. We rewrite the 
formula (12T1) for the probability amplitude ilj{m,t) into the form 

^ lik 

^{m,t) = J2 ^e^'^'^^^'^-'^' {v,{k),^)v,{k), (27) 



where we have introduced a = y. Due to the fact that we now concentrate on the 
amplitudes at the positions m ~ t we have to consider modified phases 

Uj{k) = ujj{k) - ak. (28) 

The peak occurs at such a position mo where both the first and the second derivatives of 
u}j{k) vanishes. The velocity of the peak is thus ao = Hence, solving the equations 

-v., r-1 , Vp(r + 1) cos (fc^) 

LJ-, (k) = 1 , — a = 0, 

2 ^4 + 2p [cos(A;(r + !))-!] 



uj'^ik) 



1 ^(r + 1) cos (fc^) 



v/4 + 2p [cos(A;(r + 1)) - 1] 



a 



(P-I)v/P(r + I)^sin(fc^) 
V2[2-p + pcos(A;(r + l))]2 

(29) 

for a determines the velocities of the left and right peak vl,r- The third equation is 
independent of a and we easily find the solution 

Ann , 27r(2n + 1) , , 

r + 1 r + 1 

Inserting this ko into the first two equations we find the velocities of the left and right 
peak 

r — 1 r + 1 

VL = —^ 

r — 1 r + 1 , . 

Vr = + -^VP- (31) 
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We illustrate this result in Figure [3] where we show the probability distribution generated 
by the quantum walk for the particular choice of the parameters r = 3, P = The 
initial state was chosen according to a = ^ and = vr. Since the velocity of the left 
peak vl is negative this biased quantum walk is recurrent. 



"+0 



0.08 



0.04 
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Figure 3. Velocities of the left and right peak of the probability distribution generated 
by the biased quantum walk on a line and the recurrence. We have chosen the 
parameters r = ?>, a = p = ^ and (f = n. The left peak propagates with the velocity 
Vl ~ —0.68, the velocity of the right peak is vr « 2.68. In between the two peaks the 
probability distribution behaves like while outside the decay is exponential. Since 
the velocity vl is negative the origin of the walk remains in between the left and right 
peak. Consequently, this quantum walk is recurrent. 



The peak velocities have two contributions. One is identical and independent of p, 
the second is a product of r and p and differs in sign for the two velocities. The obtained 
results indicate that biasing the walk by having the size of the step to the right equal 
to r results in dragging the whole probability distribution towards the direction of the 
larger step. This is manifested by the term ^ which appears in both velocities vl^r 
with the same sign. On the other hand the parameter of the coin p does not bias the 
walk. As we can see from the second terms entering the velocities it rather influences 
the rate at which the walk spreads. 

As we have discussed above the biased quantum walk on a line is recurrent if and 
only a Vl is negative and vji is positive. The form of the velocities (I3T1) implies that 
this condition is satisfied if and only if the criterion (l26l) is fulfilled. 

6. Mean value of the biased quantum walk and recurrence 



As we discuss in the Appendix A the classical random walks are recurrent if and only 
if the mean value of the position vanishes. We now show that this is not true for biased 
quantum walks, i.e. there exist biased quantum walks on a line which are recurrent but 
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cannot produce probability distribution with zero mean value. This is another unique 
feature of quantum walks compared to the classical ones. 

In the Appendix B we derive the following formula for the position mean value 

'X' 

1 



(l_^l_^)(a(r + l)-l) + 

Va(l-a)(l - v^r^)(l -p){r + l)cosv? 



+ 



Vp(i-p) 



r — 1 



+ ^^l-p + 0{t-'). (32) 

We see that for quantum walks the mean value is affected by both the fundamental walk 
parameters through r and p and the initial state parameters a and (p. The mean value 
is typically non- vanishing even for unbiased quantum walks ( with r = 1 ). However, 
one easily finds that the initial state with the parameters a = 1/2 and 93 = 7r/2 
results in a symmetric probability distribution with zero mean independent of the coin 
parameter p. Indeed, the quantum walks with r = 1, i.e. with equal steps to the right 
and left, do not intrinsically distinguish left from right. On the other hand the quantum 
walks with r > 1 treat the left and right direction in a different way. Nevertheless, one 
can always find for a given r a coin parameter po such that for all p > po the quantum 
walk can produce a probability distribution with zero mean value. This is impossible 
for quantum walks with p < po and we will call such quantum walks genuine biased. 

Let us now determine the minimal value of p for a given r for which mean value 
vanishes. We first find the parameters of the initial state a and ip which minimizes 
the mean value. Clearly the term on the second line in ( l32l) reaches the minimal value 
for (fQ = IT. Differentiating the resulting expression with respect to a and setting the 
derivative equal to zero gives us the condition 



(2a-l)Jp(l- p) , , 

2 + ^ = (33) 

p^y a(l — a) 

on the minimal mean value with respect to a. This relation is satisfied for = |(1 — ^/p). 
The resulting formula for the mean value reads 

x\ _r-l ^ (l-v/W-p)(l + r) ^^^^ 



1 1 aQ,ipo 2 2a/ (1 — p)p 

This expression vanishes for 

Mr)=(^)'. (35) 

Since flMl) is a decreasing function of p the mean value is always positive for p < po 
independent of the choice of the initial state. For p > po one can achieve zero mean 
value for different combination of the parameters a and ip. 

The formula (!35l) is reminiscent of the condition (!26|) for the biased quantum walk 
on a line to be recurrent. However, r is in ( l35l) replaced by r^. Therefore we find the 
inequahty p/j < po. Hence, the quantum walks with the coin parameter Pr < p < po 
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are recurrent but cannot produce a probability distribution with zero mean value. We 
conclude that there are genuine biased quantum walks which are recurrent in contrast 
to situations found for classical walks. 

7. Conclusions 

We have analyzed one dimensional biased quantum walks. Classically, the bias leading 
to a non-zero mean value of the particle's position can be introduced in two ways — 
unequal step lengths or unfair coin. In contrast, for quantum walks on a line the initial 
state can introduce bias for any coin. On the other hand, for symmetric initial state 
modifying only the unitary coin operator while keeping the equal step lengths will not 
introduce bias. Finally, the bias due to unequal step lengths may be compensated for 
by the choice of the coin operator for some initial conditions. For this reason we have 
introduced the concept of the genuinely biased quantum walk for which there does not 
exists any initial state leading to vanishing mean value of the position. 

We have determined the conditions under which one dimensional biased quantum 
walks are recurrent. This together with the condition of being genuinely biased give rise 
to three different regions in the parameter space which we depict as a "phase diagram" 
in Figure m 




r 

Figure 4. "Phase diagram" of biased quantum walks on a line. The horizontal axis 
represents the length of the step to the right r and the vertical axis shows the coin 
parameter p. The dotted line corresponds to the recurrence criterion (|26p. while the 
squares represent the condition (|35p on the zero mean value of the particle's position. 
The quantum walks in the white area are transient and genuine biased. In between 
the two curves (light gray area) we find quantum walks which are recurrent but still 
genuine biased. The quantum walks in the dark gray area are recurrent and for a 
particular choice of the initial state they can produce probability distribution with 
vanishing mean value. 

The presented results allow for generalization to biased quantum walks in higher 
dimensions assuming we keep the coin operator in a tensorial form. For non-factorizable 
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coin operators in higher dimensions it remains an open question when they are recurrent 
or transient. 
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Appendix A. Recurrence of classical biased random walk on a line 

Classical random walks on a line can be biased in two ways - the step in one direction 
is greater than in the other one and the probability of the step to the right is different 
from the probability of the step to the left (see Figure lAlj) . 



1 -p P 




Figure Al. Schematics of tlie biased random wallc on a line. Tlie particle can move 
to tlie right by a distance r with the probability p. The length of the step to the left 
is unity and the probability of this step is 1 — p. 



Consider a random walk on a line such that the particle can make a jump of length 
r to the right with probability p or make a unit size step to the left with probability 
1 — p. The random walk is recurrent if and only if the probability to find the particle 
at the origin at any time instant t does not decays faster than This probability is 
easily found to be expressed by the binomial expression 



Poit) 



tr t 

Pj r-\-lpr-\ 



t 

tr 
r+1 



With the help of the Stirling's formula 

'U 

. e 



V27rn 

we find the asymptotic behaviour of the probability at the origin 



r + 1 



. + 1 

p\ r + lpr- ' 



ip r + 1 



The asymptotics of the probability Poit) therefore depends on the value of 

.r + 1 



[1 — p) r + lpr 



f r + 1 



(Al) 



(A.2) 



(A.3) 



(A.4) 
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Since g < 1 the probability Po(t) decays exponentially unless the inequality is saturated. 
Hence, the random walk is recurrent if and only if q equals unity. This condition is 
satisfied for 

P=VTY^ (A.5) 

i.e. the probability of the step to the right has to be inversely proportional to the length 
of the step. 

This result can be well understood from a different point of view, as we illustrate in 
Figure IA2[ The spreading of the probability distribution is diffusive, i.e. cr ~ Vt. The 
probability in the a neighborhood of the mean value (x) behaves like t~2 while outside 
this neighborhood the probability decays exponentially. Therefore for the random walk 
to be recurrent the origin must lie in this a neighborhood for all times t. However, if 
the random walk is biased the mean value of the position (x) varies linearly in time, 
thus it is a faster process than the spreading of the probability distribution. In such 
a case the origin would lie outside the a neighborhood of the mean value after a finite 
number of steps leading to the exponential asymptotic decay of the probability at the 
origin Po{t). Hence, the random walk is recurrent if and only if the mean value of the 
position equals zero. Since the individual steps are independent of each other the mean 
value after t steps is simply a t multiple of the mean value after single step, i.e. 

= t{x{l)) = t [p{r + 1) - 1] . (A.6) 

We find that the mean value equals zero if and only if the condition (1A.5I) holds. 



Appendix B. Mean value of the particle's position for a quantum walk on a 
line 

In this Appendix we find the explicit form of the position mean value of the particle. 
With the help of the weak limit theorem [36] we express the mean value after t steps in 
the form 

(t) ^ E r § ^^-^^^ {v,{k),^P)v,{k), (B.l) 
j=i -^-^ ^ 

up to the corrections of the order 0{t^^). Here Vj{k) are eigenvectors of the unitary 
propagator U{k), ujj{k) are the derivatives of the phases of the corresponding eigenvalues 
and is the initial state expressed in (|T8l) . The derivatives of the phases are given in 
( l23l) . We express the eigenvectors in the form 

V2{k)=n2{r,k) (^^/rr-p^.^ + e^M'^-^'^y . 

(B.2) 




Figure A2. Spreading of the probability distribution versus the motion of the mean 

value of a biased classical random walk on a line. While the spreading is diffusive 
(cr ~ v^) the mean value propagates with a constant velocity {{x) ~ t). The probability 
inside the a neighborhood of the mean value (x) behaves like t~^. On the other hand, 
as we go away from the a neighbourhood the decay is exponential. Hence, if the mean 
value (x) docs not vanish the origin of the walk leaves the a neighborhood of the mean 
value. In such a case the probability at the origin decays exponentially and the walk 
is transient. 



x\^ f ^ ^ fia,^,P,r,u)du _ + (B.4) 

J 2(r + IjTT 1 + ,/pcosMi 1 - ->/psinM2 



The normalization factors of the eigenvectors read 

ni{u) — 2 — 2y^cos {u — arcsin [y^sinw]) , 
712 (li) = 2 + 2^/pcos {u + arcsin [y^sinii]) , 

(B.3) 

where we denote u — Mz±l) to shorten the notation. The mean value is thus given by 
the following integral 

{r+l)7r 

2(r + l)7r [l + y^cosMi] [l - ^/psmu2] 

u 

where 

ui — u + arcsin(^/p sin u), U2 — u + arccos(^/p sin u) , (B.5) 
and the numerator reads 

/(a, (/?, p, r, m) = (1 - p) [r - 1 + p (a + r{a - 1)) (1 + cos(2-u)) + 

+ ^/a{l — a)^/p{l — p){r + 1) (cos + cos((^ + 2u)) . 

(B.6) 

Performing the integrations we arrive at the result 

(l-yW)(a(r + l)-l) + 
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^ v/a(l-a)(l - v/r^)(l - p){r + l)cos<^ 



Vp(i-p) 



r — 1 



+ ^V^ + Or^). (B.7) 
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